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Abstract 

Three Large Eddy Simulations (LES) for a lean-direct injection (LDI) combustor 
are performed and compared. In addition to the cold flow simulation, the effect of 
radiation coupling with the multi-physics reactive flow is analyzed. The flamelet 
progress variable approach is used as a subgrid combustion model combined with 
a stochastic subgrid model for spray atomization and an optically thin radiation 
model. For accurate chemistry modeling, a detailed Jet-A surrogate mechanism is 
utilized. To achieve realistic inflow, a simple recycling technique is performed at the 
inflow section upstream of the swirler. Good comparison is shown with the experi- 
mental data mean and root mean square profiles. The effect of combustion is found 
to change the shape and size of the central recirculation zone . Radiation is found 
to change the spray dynamics and atomization by changing the heat release dis- 
tribution and the local temperature values impacting the evaporation process. The 
simulation with radiation modeling shows wider range of droplet size distribution by 
altering the evaporation rate. The current study proves the importance of radiation 

modeling for accurate prediction in realistic spray combustion configurations, even 
for low pressure systems. 
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1 Introduction 


All modern combustors must balance the need for stability and performance 
with the goals of efficiency and emissions reduction. The emissions from sta- 
tionary as well as non-stationary gas turbines have a profound effect on the 
environment and the balance of the eco-system [1], For example, green house 
gases like Carbon dioxide C0 2 , water vapor H 2 0, and methane CH 4 are pri- 
mary by-products from the burning of fossil fuels in internal combustion and 
gas turbine engines. These gases contribute to the global warming effect by 
absorbing the sun infrared rays. Other important pollutants are CO and NO x . 
While CO has known hazardous health effects, NO causes acid rain and de- 
pletes the ultra-violet protective layer known as the ozone layer [2], Finally, 
hydrocarbon fuels under rich conditions generate soot, which is known for its 
carcinogenic precursors and associated organics [3,4], A common technique to 
reduce emissions of NO^, CO, soot, and unburned hydrocarbons is to operate 
at lean premixed conditions, which in return leads to stability problems [5- 
7]. For aircraft engines, this poses a safety concern and is therefore typically 
avoided. Modern engines for next generation supersonic transportation air- 
craft, however, face even more stringent emissions regulations. Because of 
the high cruise altitude of these aircrafts in the lower stratosphere, engine 
emissions contribute directly to depleting the ozone layer [8]. To meet this 
challenge, several techniques have been proposed that are based on essentially 
premixed combustion. Among these are lean- premixed-prevaporized combus- 
tion (LPP) and lean direct injection combustion (LDI [9]). The latter is the 
focus of the present study. 
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The numerical prediction of the emissions in each system is faced by many 
challenges [10,11]. Gas turbine engines operate under turbulent conditions us- 
ing complex fossil fuels [12] that introduce a wide range of length and time 
scales [13]. The resolution of these scales usually require fine mesh [14], which 
can lead to costly computation. In addition, the production of pollutants usu- 
ally involves a large number of chemical reactions that expands the numerical 
effort. Furthermore, the complexity of the system is exacerbated as most of 
gas turbine engines fuels are injected in the liquid phase. The breakup of the 
spray droplets and their subsequent evaporation is another complex process 
that is not well understood yet [15,16] and which might be inherently coupled 
with the formation of pollutant and other physical process such as radia- 
tion [17]. Therefore, efficient numerical tools are needed to simulate reactive 
multiphase flow in complex geometries under high Reynolds number turbulent 
conditions [18-20]. 

In the past decade, large eddy simulation (LES), has proved to be a promising 
tool for turbulent reactive flow simulations [20] using both structured [21-26] 
and unstructured meshes [27-30]. In the current work, the full single element 
Lean Direct Injection (LDI) combustor is simulated in the LES framework. 
The Lean Direct Injection configuration is a good candidate for low emissions 
gas turbine requirements, where the combined geometry of a swirler and a 
Venturi nozzle results in stable combustion with potentially ultra-low NO x and 
soot emissions. Air is injected in a swirler and it mixes with the fuel droplets to 
atomize, break up aero dynamically and partially premix in the Venturi before 
entering the combustor. The combination of the swirler and the Venturi has 
proven to maximize the atomization performance and minimize pressure drop 
across the injector [31]. The Venturi nozzle also provides sufficient residence 

time for the fuel droplets to vaporize and mix uniformly with the swirled air 
droplets in a lean mixture. 
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Here, non-reactive and reactive multi-physics simulations will be shown for an 
LDI system [32-34], The current paper presents simulations that couple the 
interaction of turbulence with the important physical processes such as spray 
dynamics, chemistry, and radiation in a realistic configuration. The combustor 
emission characteristics, radiation coupling, flow/flame interaction, as well as 
the spray characteristics will be shown and compared with experimental data. 
A description of the modeling approach is outlined in the next section. 

For accurate chemistry modeling, a detailed Jet-A surrogate mechanism is 
utilized [35]. The number of species and reactions considered are 122 and 900, 
respectively. Jet-A is known to be a common fuel for gas turbine engines. In 
the context of LES, the detailed Jet-A surrogate mechanism will be introduced 
in the frame work of the flamelet approach [36]. 


2 Combustion and radiation modeling 


The combustion model used is the flamelet /progress- variable (FPV) approach [37,38]. 
The FPV approach is based on the flamelet concept, which relates the species 
mass fractions and energy to the mixture fraction through the flamelet equa- 
tions 


<9$ x <9 2< £ 

P ' ~dt ~ P 2 dz* = u ’ 


( 1 ) 


where Z is the mixture fraction, p is the density, $ = (T, is a vector 
that contains the species mass fractions yi and the temperature T. uj is the 
source-term vector and x is the scalar dissipation rate, which represents the 
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local rate of molecular mixing and which is strongly impacted by the large 
scale turbulent mixing process. The steady laminar flamelet model (SLFM) 
tabulates the solution of Eq. (1) (without the unsteady term) as a function of 
the mixture fraction and the stoichiometric scalar dissipation rate y st . Hence 
an equation of state (EOS) relating reactive scalars to the mixture fraction 
and its dissipation rate is written as 

$ = (Z, Xst) • (2) 

The solutions to the flamelet equations are characterized by the non-linear in- 
teraction between chemistry and transport. For increasing dissipation rate, the 
flame temperature decreases, which ultimately leads to extinction. In SLFM, 
only non-extinguished solutions of the flamelet equations are considered, which 
limits the applicability of the model. It was shown by Pierce and Moin [37], 
for example, that SLFM predicts too much heat release in the essentially 
non-reactive regions close to the nozzle for a flame that is aerodynamically 
stabilized by a recirculation region. To overcome this limitation, Pierce and 
Moin [37] proposed a tabulation based on a reaction progress parameter A 
instead of the scalar dissipation rate, which leads to the Flamelet/Progress 
Variable (FPV) formulation of the steady flamelet model. 

As a result, the new state vector is given as T = (Z, A). The reaction 

progress parameter A is defined according to Ilime et al. [38] through a reaction 
progress variable C such that it is independent of the mixture fractional. 
The progress variable C can be a combination of certain species representing 
combustion products, temperature or sensible enthalpy. The progress variable 
here is the sum of the mass fractions of CO 2 , CO, H 2 0, and H 2 . The reaction 
progress parameter A is then defined for each flamelet as the stoichiometric 
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value of the reaction progress variable. Hence, the reaction progress parameter 
an be obtained from an inversion of the fla /m el et library. 

C = T c (Z, A) , (3) 

assuming that inversion is unique, which puts a restriction on the definition 
of the reaction progress variable. 

The state vector $ = (Z, A) is obtained by the solution of Eq. (1), which 

includes all possible steady state solutions. Then, a joint probability density 
function (PDF) needs to be evaluated to obtain the Favre filtered (~) values 
of the state vector. Since Z and A are defined to be independent variables, by 
using Bayes' theorem, the joint PDF at a certain time and position in space 
is given as 


P (Z, A; x, t ) = P (A|Z; x, t) P (Z; x, t) , (4) 

where marginal PDF of the mixture fraction P (Z) is assumed to be a beta 
distribution [37], and the conditional reaction progress parameter PDF is a 
delta function. Since the (3 - distribution is characterized by the scalar mean Z 
and the scalar variance Z" 2 , and the delta-function by the mean A the flamelet 
tabulation can be written in the form: 

* = As (z > ^7 d ■ (5) 

With Eq. 5 and 

C = T A (z,r \\) , ( 6 ) 
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the mean reaction progress parameter as an independent parameter in the 
tabulation can be replaced by C, which is readily available in a simulation 
through the solution of a filtered transport equation for C . As previously 
mentioned, for accurate temperature predictions, radiation might have to be 
considered. Ihme and Pitsch [39] have extended the above FPV approach to 
account for radiation as an unsteady process. In this model, the unsteady 
flamelet equations are solved considering radiative heat losses. The unsteady 
solution starts from the steady adiabatic solution and proceeds until a new 
steady state is found. The time can be replaced by the enthalpy, which is 
defined to include the heats of formation, and which therefore appears as a 
new parameter in the tabulation. Consequently, the extended FPV EOS can 
be written as 

n = ^{z,z^,c,H ) . (7) 

As a result, four parameters have to be computed as part of the LES so- 
lution, namely the hltered mixture fraction Z, the subfilter scalar variance 
Z " 2 , the hltered progress variable C, and the hltered enthalpy H. Here, the 
scalar variance Z" 2 is evaluated using an algebraic model [37], which assumes 
homogeneity and local equilibrium for the subgrid scales and is given by 

pZ 7r2 = C z A 2 p\VZ\ 2 , (8) 

where the mixture fraction coefficient Cz is computed dynamically and p is 
the hltered density. Finally, three additional transport equations are solved 
for Z, H and C as follows: 

^ + V • (p8Z) = V • (pBVZ) + V • tT* , (9) 
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( 10 ) 


^ + V- [puc) = V • (pDWC) + V ■ T^ es + PUJC , 

^ + V ■ (puH) = V ■ (pDVtf) - ~q r , (11) 


where D is the filtered molecular diffusivity, assuming equal species diffusivi- 
ties and unity Lewis numbers, and ioc is the hltered chemical source term of 
the progress variable, which is tabulated in the flamelet library. In the above 
equations, the residual stresses r res are modeled by a dynamic approach [37]. 
For an optically thin radiation model [40], the hltered radiation source term 
is 


q r = 4a (T‘ - Z‘_) ^(pia,) . 


(12) 


where pi and ai are the partial pressure and Planck mean absorption coeffi- 
cient of species i, respectively, a is the Stefan Boltzmann constant, and 
is the ambient reference temperature at one atmosphere. To account for the 
unsteadiness of NO^ production, the hltered mass fraction of NO, Yno is com- 
puted by solving a transport equation rather than extracting the values from 
the hamelet tables [39]. The hltered LES equation can then be written as 

+ V ■ (pu Fno) = V ■ (pDVY no) + pch no , (13) 


where the hltered production rate Ono is closed according to the model by 
Ihme and Pitsch [39]. More details about the FPV and the radiation model 
can be found in [39]. 


NASA/TM— 2011-217111 


3 Secondary breakup spray model 


Jet-A fuel spray droplets are injected from the centerline simplex atomizer, 
with a specified drop-size distribution. The hollow cone spray has a total angle 
of 90°. The spray models employed have been described in detail by Apte et 
al. [41,28,42], Here, we summarize the most important features. 

A stochastic subgrid model for spray atomization is used here. The initial 
liquid film is approximated by a droplet of size equal to the nozzle diameter. 
The model views the subsequent particle breakup of the injection of a spec- 
ified drop-size as a discrete random process, where droplet size is treated as 
a time-dependent stochastic variable. The temporal and spatial evolution of 
the particle size distribution is then governed by the Fokker-Plank differential 
equation. This evolution leads to the formation of new droplets and destruc- 
tion of the parent ones. The new droplets position and velocity vectors are 
tracked by a Lagrangian algorithm in the physical space. Two-way coupling 
between the gas phase and the liquid droplets is achieved by the source terms 
in the gas-phase equations, which represent the mutual effect of mixing and 
momentum/energy transport. These effects are induced by droplet breakup, 
evaporation, and convection by the flow-field. Simultaneously with the process 
of secondary droplet breakup and mixing with the gas phase, the liquid phase 
evaporates under the appropriate conditions. The evaporation model solves 
a set of ordinary differential equations (ODEs) that describe the variation 
of the droplet mass and temperature due to evaporation and assumes that 
non-equilibrium effects inside the droplet volume are neglected. A third-order 
Runge-Kutta scheme is used to integrate the set of ODEs with the minimum 
physical drop time scale. The mathematical formulation of the above scheme 
is omitted here for brevity and can be found elsewhere [42], 
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4 Numerical setup 


The simulation is performed by the unstructured LES code CDP. CDP is 
a set of massively parallel unstructured finite-volume flow solvers developed 
specifically for large eddy simulation. The solver used to perform the reactive 
simulations reported in this work is a node-based low-Mach number solver [43] . 
A the fully incompressible solver is used for the non-reactive case for compu- 
tational efficiency. The LES governing equations and the sub-grid momentum 
closures are omitted here for brevity. 

The experimental setup and data are provided by Farhad et al. [44], The ge- 
ometry of the single-element combustor is shown in Fig. la with the inflow 
pipe upstream of the swirler. Liquid fuel is injected through the center (closed 
for the non-reactive case), while air is injected through a swirler with vane an- 
gles of 60°. The global swirl number is 1. The swirler is composed of six helical 
vanes with an effective area of 870 mm 2 . The fuel droplets from the centerline 
and the swirled air are mixed in a converging-diverging Venturi nozzle. The 
helical air swirler has an inside diameter of 9.3 mm and an outside diameter of 
22.1mm. The combustor has a square cross-sectional area of 50.8 x 50.8 mm. 

The mesh uses 1,072,640 hexahedral elements, which is slightly above the 
861, 823 element mesh used in the RANS calculations [44,45] and about half 
the LES structured mesh used by Patel and Menon [22], The unstructured 
mesh is shown in Fig. 1. The mesh distribution along the z = 0 plane and 
in the vanes is shown in Fig. lb, while the y — z plane grid distribution 
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is presented in Fig. lc. The mesh cells are clustered toward the centerline 
and toward the walls of the combustor. Here, we extend the original RANS 
mesh [44] to account for more realistic inflow. 

An inflow bulk velocity of 20.14 m/s is provided through a tube upstream of 
the swirl injector. The inflow air is at temperature T a = 294 K and pressure 
of 1 atm, while the global equivalence ratio is 0.75. The fuel spray has an 
inflow mass flow rate of 4.15.10—04 kg/s and an initial velocity of 20 m/s. 
The droplets are injected with a Rosin-Rammler distribution 


where F m (D) is the cumulative distribution function of the droplet diameter 
.D, and the parameters n = 1.34 and 6 = 24 /i m are chosen to curve-fit 


turbulence. The recycling is done until the inflow profile of the upstream pipe 
recovers a realistic inflow turbulent boundary layer. The boundary conditions 
used are inflow/outflow in the ^-direction and adiabatic walls in the y and z 
directions. 

Statistics are collected over a physical time of 0.032 seconds after the initial 
transient for the reactive flow cases and over 0.167 seconds for the cold-flow 
case. Based on the bulk inflow velocity and the full length of the combustor, 
this represents two complete flow-through times for the reactive cases and 
about seven flow-through times for the non-reactive case. The simulation was 
performed on a Linux cluster at Stanford University with inhniband intercon- 
nection and dual Intel Clovertown (Quad Core) 2.33 GHz processors with 8 GB 



(14) 


the distribution of the experimental data at the first measured location (x = 
3mm). An inflow recycling technique is used [46] to achieve realistic inflow 
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RAM per node. With the current mesh resolution the code scaling was found 
to level off at 96 processors for the spray-combustion case. The computational 
time per flow through time was 900 CPU hours for the reactive case. 


5 Results and discussion 

In this section, two LES simulations are compared to the experimental data. 
The statistics for the two simulations, with and without radiation, are collected 
over the same total run time. Unless mentioned otherwise, only the case with 
radiation is considered for the comparison with the radial experimental pro- 
Eles. Almost all available experimental data are used in the comparison and 
selected non-reactive flow data are included to show the effect of heat release. 

First, the time-averaged flow features of the reactive flow will be shown and 
compared with the non-reactive case with an emphasis on the radiation effect. 
Subsequently, the mean and root mean squared (nns) values of the gas flow- 
field, the scalars, and the spray characteristics will be compared with the 
experimental data. 

Due to the swirling flow, an adverse pressure gradient is generated in the 
axial direction and a vortex breakdown bubble (VBB) is established. The 
LDI is characterized by a central recirculation zone (CRZ) in addition to 
recirculation zones (RZs) at the corner of the dump plane. Three RZs had 
been observed previously in the non- reactive case [46,47,45]. Figures 2a and 2b 
show isosurfaces of mean axial velocity for ( U ) = Orn/s for both the reactive 
and non-reactive cases, where the brackets denote a time average. The length 
of the CRZ in the non-reactive case is approximately twice as that in the 
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reactive one. However, the reactive CRZ has a larger effective bulk diameter. 
Also, a small toroidal RZ close to the divergent part of the Venturi exists in 
the reactive flow case as shown by the side view in Fig. 3, which is caused 
by flow separation. A more compact and stronger reactive CRZ is due to the 
expansion of the flowheld by heat release effects as shown later. 

Figure 4 shows vector plots for the time-averaged total velocity magnitude for 
the reactive flow case with radition. Only the projections in the y = 0 plane 
are shown here for clarity. Figure 4a shows that the CRZ starts early inside 
the divergent part of the Venturi nozzle, in addition to a small RZ observed 
just downstream of the injector exit. These two small RZs inside the Venturi 
enhance mixing and atomization of the spray with the swirling air. However, at 
the injector exit, the droplets and the spray sheet have enough momentum to 
penetrate this RZ toward the Venturi exit as shown in Fig. 4b. This observation 
will be confirmed later by the axial and radial velocity profiles. Figure 4c shows 
the velocity magnitude vector plot downstream of the dump. The CRZ has 
an azimuthal vortex at the location x — 17 mm downstream of the dump 
plane. Figure 4c also shows an additional annular vortex at the burner corners 
where the second RZ is located. Corner vortices are found to extend along the 
combustor wall corners up to the exit plane in the y — z plane (not shown 
here), which will affect the spray distribution as shown later. 

An instantaneous temperature distribution is shown in Fig. 5. Figure 5a shows 
a snapshot of the y = 0 plane. Figures 5b-d show three consecutive snapshots 
at the location x = 11 mm. Figures 5a and 5d are at the same instant in time. 
The temperature shows high gradients at the shear layer location and around 
the CRZ, which leads to strong mixing. At the edge of the CRZ, fresh cold 
gases (dark areas) are engulfed into the CRZ. Inside the core of the RZ, the 
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temperature shows nearly uniform distribution due to the strong mixing the 
fluid undergoes before entering the CRZ. 

A comparison of the heat release distribution is shown in Fig. 6. The flame 
without radiation shows locally lower mean values than the case with radia- 
tion, for which the heat release is more distributed. The reason is that radiation 
causes the temperatures to be locally lower, which leads to smaller evapora- 
tion rates of the liquid fuel and consequently a lower fuel vapor mass fraction 
as observed in Fig. 7. 

Two snapshots of the spray distribution from the two simulations are shown 
in Fig. 8. The spray exhibits a 90° cone angle. Some of the spray droplets 
are trapped inside the corner vortex, which results in an increase in the mean 
droplet diameter close to the walls. The simulation with radiation shows a 
wider range of droplet sizes. For the shown instantaneous snapshots with ra- 
diation, droplet size ranges between 0.5 /mi and 50 /mi, while without radiation 
the range is between 0.7 /mi and 40 /mi. This is attributed to a higher evapo- 
ration rate when radiation is not included. Since no primary breakup model is 
employed here, the spray spectrum shows an intermediate regime, where the 
spray injected from the nozzle undergoes secondary breakup, followed by a di- 
lute regime where the droplets evaporate. The intermediate regime is mainly 
inside the Venturi, while the dilute regime extends until all the droplets are 
evaporated. 

The steady state droplet and parcel-size histograms for the whole combustor 
volume for both cases are shown in Fig. 9. In general, the case without radi- 
ation has a lower number density of droplets over the whole spectrum. This 
will subsequently affect the flow dynamics and the flame structure (as shown 
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earlier). The radiation effect on droplet evaporation can be assessed in Fig. 10. 
This figure shows the time-averaged rate of droplet mass evaporation at the 
plane y — 0. This value represents the time average of the source term of the 
droplet mass variation. The case without radiation shows higher evaporation 
rates in the vicinity of the injector and in the flame surface surrounding area 
due to the reasons discussed earlier. 

The LES simulation with radiation modeling is compared with the experimen- 
tal data next. Figures 11a and lib show spray droplet mean diameter Dio and 
Sauter mean diameter D 32 . Overall good agreement is observed. At the first 
three axial locations a peak value is shown that starts at a radial distance of 
8 mm from the centerline (for x = 3 mm) and then spreads radially towards 
the wall (for x > 3 mm). This location is an indication of the spray sheet 
cone angle. Since the spray droplets directly emitted from the injector have 
a large droplet size, the profiles show peak values at these radial locations. 
Downstream of x — 15 mm some droplets are trapped in the wall-corner vor- 
tex. It is interesting to note that the value of the mean diameter is increasing 
downstream due to evaporation of the smaller diameter droplets (larger than 
the cut-off value of the experimental measurements). This is consistent with 
the experimental data. 

The spray droplets' mean axial velocity distribution is shown in Fig. 12. Good 
agreement is observed for all locations. Consistent with Fig. 11, the LES shows 
over-predicted droplet mean axial velocity at the locations of over-prediction 
in the droplet size. Again, over-prediction is expected at the first few locations, 
where the small-size droplets are not considered in the experiment. 
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To show the extent of the centerline RZ and the effect of heat release and 


radiation on the flow features, the centerline mean axial velocity profiles are 
shown in Fig. 13. Reactive and non- reactive data from experiments are com- 
pared with the LES results with and without radiation. The figure shows that 
the non-reactive RZ extends to about 100 mm downstream of the dump, while 
the reactive RZ is stronger and more compact (about 50 mm in length) . The 
expansion of the flow by heat release changes the local velocity and pressure 
distribution. This local expansion (in the flame region) increases the pressure 
gradient and leads to faster flow deceleration. This results in a shorter but 
wider RZ. The CRZ without radiation is over-predicted by about 5 mm due 
to the change in the flame structure and the spray distribution as shown in 
Fig. 6. 

As discussed below, the LES simulation shows that inside the divergent section 
of the Venturi a small RZ is established, followed by an increase in the velocity 
due to a sudden expansion after the dump. This is followed by a steep reduction 
inside the RZ due to the adverse pressure gradient. There are no experimental 
data inside the Venturi, but the prohle captures the initial peak after the dump 
accurately, which confirms the foregoing physical conclusions. Note that the 
current simulation captures the centerline prohle initial peak. The deviations 
near the injector will be discussed later. 

The mean and rms axial velocities are shown in Figs. 14 and 15. The first 
three locations (x = 3 mm, x = 5 mm and x = 9 mm) are found to be the 
most difficult to match with the experiment. These locations exhibit high un- 
steadiness and high measurement errors. The experimental measurements [34] 
use a Phase Doppler Particle Analyzer (PDPA). A problem usually related to 
the measurements is how to distinguish between the seed particles to measure 
the gas phase velocity and the spray droplets. In the LDI experimental mea- 
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surements [34] the nominal diameter of the seed particles used was 1 /irn. To 
distinguish between the spray and the seed particles, all particles smaller than 
4 /irn were used to represent the gas-phase velocity. However, inside the CRZ 
the small spray droplets just ejected from the nozzle still have high momen- 
tum in the positive direction. As a result, the experimental velocity profiles 
near the exit (i.e., x = 3 mm, x= 5 mm and x = 9 mm) are all positive despite 
the existence of the CRZ. Therefore, at these locations we expect to under- 
predict the experimental data. A previous RANS simulation [45] shows also 
higher deviation at x = 3 mm, and comparable results at x = 9 mm. The LES 
simulation by Patel and Menon [22] also shows lower velocity at the location 
x = 5 mm. Starting from the x — 15 mm location, the current LES simu- 
lation shows good results compared with the experimental data. The CRZ 
extends to approximately 50 mm downstream of the dump. At the x = 92 mm 
location, the prohle shows a nearly uniform velocity that recovers the inflow 
mass flow rate at the exit cross-sectional area. This uniformity is important 
in gas turbine engines to minimize the thermal stresses on the turbine blades. 
A slight deviation is found towards the walls. This can be attributed to the 
fact that the CRZ in the current simulation is wider than in the experiment. 
As a result, the flow is more compressed between the CRZ and the walls and 
the velocity is higher. 

The rms axial velocity prohle is shown in Fig. 16. The rms peaks close to 
the dump plane at the onset of the shear layer and then decays axially as 
the shear layer decays downstream. The predicted and experimental axial rms 
profiles show good agreement. The high rms values at the first three locations 
confirm the high unsteadiness of the shear layer. For example, at x = 5 mm 
the rms value is around 100% of the corresponding mean value. At the first 
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three locations, over-prediction is noticed at the y — 15 mm radial location. 
The discrepancy with the experimental data at the first three locations is due 
to the reasons mentioned in the previous paragraph. The LES data predict 
two rms peaks that start at the radial location y = 9 mm at the first axial 
location. These two peaks spread downstream with the spread of the shear 
layer. Downstream of x — 20 mm, the peaks start to move again toward the 
centerline following the CRZ surface. After reaching the x = 46 mm location, 
the peaks are completely merged. After this location, no reverse flow exists 
and the shear layer is completely diffused. Similar observations apply to Fig. 
17, where the radial profiles of the turbulent shear stress are plotted. The 
shear stress shows antisymmetric radial profiles. Good agrement is shown, with 
deviations at the first two locations. The zero shear stress at the centerline 
indicates flow symmetry around the centerline. 

Figure 18 shows the mean velocity component in the ^-direction at different 
axial locations. The flow is rotating in the clockwise direction (CCW) viewed 
from the outflow cross-section. The global swirl number is of the order of 1. In 
contrast to other velocity components, the results show good agreement with 
the first location at x = 3 mm. The flow shows solid body rotation around the 
centerline downstream of the location x = 29 mm. Upstream of this location 
the flow is rotating around the central RZ. 

The focus of the following discussion are the LDI emission characteristics. For 
supersonic vehicles, the NO emissions directly deplete the stratospheric ozone 
layer. Therefore, accurate models are required to predict gas turbine emissions. 
As discussed earlier, NO production is a slow process and it is coupled with 
other processes that have long characteristic time scales such as radiation. 
Radiation affects NO in an indirect way through the change in temperature, 
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which then changes the production rate. As a result, the first step is to predict 
temperature accurately. Figure 19 shows the centerline mean temperature with 
and without radiation modeling. This figure shows that without radiation the 
temperature is over-predicted. This is consistent with earlier observations. 

The radiation model employed here is an optically thin model [39]. The model 
assumes that the medium is non-absorbing and neglects the heat lost or gained 
by wall absorption and reflection. Figure 20 shows radial temperature profiles 
at the locations where experimental data are available. The temperature is 
under-predicted at the centerline. This may suggest a need for a more accu- 
rate radiation model that accounts for the absorption within the flame, the 
scattering by the media (although might be negligible in this case), and the 
wall effect. 

Figure 21 shows the centerline species distribution for CO and NO. The CO 
level along the centerline is captured reasonably well. An initial peak is ob- 
served in the rich combustion region followed by a decay in the CO level as it 
is oxidized to produce CO 2 . The centerline NO profile also shows good results. 
The effect of radiation generally improves the results. 


6 Conclusions 

Large eddy simulations for a liquid fuel injection LD1 are performed. The 
simulations show that LDI is a good candidate for low emissions, yet stable, 
gas turbine engines. Radiation modeling is found to be essential for accurate 
NO and temperature predictions. By changing the heat release distribution, 
radiation alters the spray dynamics and the droplets size distribution. Thus, 
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radiation can be important even for low pressure systems where radiation 
from flame region is relatively small, but effect of radiation on spray dynamics 
significantly affects temperature distribution. The simulations show the cur- 
rent ability of LES as a numerical tool to simulate multi-physics problems in 
complex geometries using realistic chemistry models. 
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(c) y-z mesh plane 

Fig. 1. The LDI single element geometry and unstructured hexahedral mesh. 
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(a) Reactive 


(b) Non-reactive 


Fig. 2. Iso-surface of the time-averaged axial velocity (( U ) 


0 m/s). 
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Fig. 3. Side view (y = 0 plane) to the CRZ for the reactive flow case. 
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(a) Zoom of the injector area 


(b) Zoom of droplet distribution in the 
injector area 



(c) The vortices around the CRZ and at 
the corners. 

Fig. 4. Vector plots of the mean velocity magnitude projected in the y = 0 plane 
for the reactive flow case with radiation. 


NASA/TM— 2011-217111 


24 



(c) z-y plane (d) z-y plane 

Fig. 5. Instantaneous temperature distribution for the case with radiation, (a) for 
the z-x plane, (b), (c) and (d) for the z-y plane at three consecutive time snapshots 
at x = 11 mm. 
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(a) With radiation (b) No radiation 

Fig. 6. Distribution of the time averaged heat release ( Qht )■ 



(a) With radiation (b) No radiation 

Fig. 7. Cross sectional view through the swirler at y = 0 plane for distribution 
of the time-averaged gas-phase fuel mass fraction (Yp) combined with the RMS 
temperature contours ( T rms ). 
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(a) With radiation (b) No radiation 

Fig. 8. Instantaneous droplet distribution sized by droplet mass combined with the 
instantaneous temperature distribution at the y = 0 plane. 
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Fig. 9. Radiation effect on the droplet parcel histogram. (With radiation 
without radiation — ) 
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Evap.Rate 



(a) With radiation (b) No radiation 

Fig. 10. Time-averaged droplet mass evaporation rate in (kg/s) at the central y = 0 
plane. 
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Fig. 11. Radial profiles of droplet arithmetic mean and Sauter mean diameters 


(Experiment data o, computation — ) 
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Radial distance (mm) 



(a) x = 3 rrirn-x = 15 mm 


(b) x = 20 mm-x = 46 mm 


Fig. 12. Radial profiles of mean droplet axial velocity (Experiment data o, compu- 


tation — ) 
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Fig. 13. Comparison of time-averaged centerline axial velocity (U). Non-reactive 
flow (Experiment □, computation — ), reactive flow (Experiment o, computation 
with radiation — • — , computation with no radiation • ■ ■). 
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(a) x = 3 mm - x = 20 mm (b) x = 29 mm - x = 92 mm 

Fig. 14. Mean axial velocity (U) at different axial locations (Experiment data o, 
computation — ). 
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(a) x = 3 rrirn-x = 20 mm 


(b) x = 29mm-x = 92 mm 


Fig. 15. RMS axial velocity U ims at different axial locations (Experiment data o, 
computation — ). 


NASA/TM- 


- 2011-217111 


34 


20 - 


r. 

i\ 


3 10 7o 


>° \ 


' ^ooooo 1 


/A 


/ \ 


Ooo0o °oooooo66od 


20 40 60 

Axial distance (mm) 


80 100 


Fig. 16. RMS axial velocity along the centerline (Experiment data o, computation — ). 
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Fig. 17. Radial profiles of the shear stress uv 


NASA/TM— 2011-217111 


36 



1 1 1 1 1 ' 1 ' I 

JUUJUI. 11 ii^nnnnnnnr>g*« ; 9eCt5&0CXX3000^\ J 

I 1 1 

1 i 1 i 1 i 1 i I 


20 


1 ' X=4mm ' 1 

1 I ^ 

ooooeoCSoo \ 

7 

20 

i I i 

i 1 i 1 

I 1 i_ 

- 


20 

0 

-20 

20 

0 

- 20 - 


^ l T I T l T l T l ^ 




20 F - 
0 L 
-20 L 

20 F - 
0 1 
-20 L_ 
20 
0 - 
-20 L_ 


_ X = 29mm , » » nr 

V — £^°®douo<W6ooooooou'^ 0 ^ UU 


atsniWrnnQect 


^22Q ooonu c 


X=38mm' 




irV^nn^oOtoaot 


v g22 fl a oopuu i 


, oOOUJ(fe)OoO ; 

1 I n_ 
X=46 


l ® 00 ® cr uooN 


^uuoyocbOu' 



20 
0 - 

-20 ^_ 
20 
0 - 

-20 ^_ 

-30 


X=72mm fi 
^Onnnnnnnnru.ir -nnQ* 


X=92mm. 

^^^T i nnnrnnnnr i r i r " r '' rinrv ^ or 
_i I i I i L 


SUeceSooS 0 ^oooooowo 0 


20 


Radial Distance (mm) 


Radial distance (mm) 


(a) x = 3 rrirn-x = 20 mm (b) x = 29 mm- a; = 92 mm 

Fig. 18. Mean velocity component in the y direction ( Uy ) (Experiment data o, 
computation — ) 
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Fig. 19. Comparison of the centerline mean temperature (T). (Experiment o, com- 
putation without radiation (— ) and with radiation (• — •). 
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Fig. 20. Mean radial temperature profiles ( T ) (Experiment o, computation — ) 
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Axial distance (mm) 


Fig. 21. Time averaged centerline species mass fractions profiles in ppm ((Tco) and 
(V no)) without radiation (— ) and with radiation (• — ■). 
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